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Abstract 

Linearized stability of incompressible viscous fluid flows in a thin spherical shell is studied 
by using the two-dimensional Navier-Stokes equations on a sphere. The stationary flow on the 
sphere has two singularities (a sink and a source) at the North and South poles of the sphere. 
We prove analytically for the linearized Navier-Stokes equations that the stationary flow is 
asymptotically stable. When the spherical layer is truncated between two symmetrical rings, 
we study eigenvalues of the linearized equations numerically by using power series solutions 
and show that the stationary flow remains asymptotically stable for all Reynolds numbers. 

1 Introduction 

The Navier-Stokes (NS) equations for an incompressible viscous fluid are the fundamental governing 
equations of fluid mechanics. In many cases, exact solutions can be constructed to these equations 
jni and spectral and nonhnear stability of these exact solutions can be analyzed |Hj. Our work 
addresses stability of exact solutions for the NS equations in spherical coordinates. 

The three-dimensional NS equations in a thin rotating spherical shell describe large-scale atmo- 
spheric dynamics that plays an important role in the global climate control and weather prediction 
|171 118j (see also review in [12 ). It was rigorously proved by Temam &: Ziane j22| that the av- 
erage of the longitudinal velocity in the radial direction converges to the strong solution of the 
two-dimensional NS equation on a sphere as the thickness of the spherical shell goes to zero. The 
latter model has been used in geophysical fluid dynamics since middle of the last century |16j . 

The treatment of the geometric singularity in spherical coordinates has for many years been a 
difficulty in the development of numerical simulations for oceanic and atmospheric flows around the 
Earth. Blinova jSl 1^ represented solutions in the inviscous case by the eigenfunction expansions in 
spherical harmonics. Vorticity equations were considered by Ben-Yu with the spectral method 0. 
More recent work of Furnier et al. ^Tj applied the spectral-element method to the axis-symmetric 
solutions (see |14| I19[ l23] for other applications of the spectral methods in spherical coordinates). 
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Finally, point vortex motion on a sphere was modeled by ordinary differential equations for vortex 
centers in Boatto & Cabral [5^ and Crowdy [7]. 

We address the three-dimensional NS equations for an incompressible viscous fluid, 

V-u = 0, xen, teR+, (1.1) 

u\t=o = uo, X e n, 

in a thin spherical shell ^7 = {x G M'^ : 1 < |x| < 1 + e} with e — > 0, subject to the boundary 
conditions 

u-n = 0, Vuxn = 0, x G 90. (1.2) 

Here u : $7 x is the velocity vector, p : i^x M^, M is the ratio of the pressure to constant 

density, z/ is the kinematic viscosity, n is the normal vector to the boundary dO, of the spherical 
shell Q and uq : is a given initial condition. Although Coriolis and gravity forces may be 

dynamically significant in oceanographic applications, our model is considered in a non-rotating 
reference frame and without external forces. The effects of rotation and gravity can be included 
into the model but they do not substantially alter the physical picture that emerges from the NS 
equations 

We employ the spherical coordinates {r,9,(j)) with the velocity vector u = UrSr + ugeg + w^e^, 
where (e^, eg, e^) are basic orthonormal vectors along the spherical coordinates. For completeness, 
we reproduce the three-dimensional NS equations in spherical coordinates ^J: 

du^ dur U0 dur Ud> dur ^6 ^ ^1 dp ( . 2 dur 2u 

^Ur— \ tSt^T^ - = -—^V\l\Ur^ h 



dt dr r 89 rsinO dcp r dr \ r dr 

dug due Uq dug U^ dug Uj-Ug cot 9 I Qp 

dt ^ dr r d9 r sin 9 dip r r r d9 

2 dur ug 2 cos 9 duj 



r 



+ 



du^ ^ du^ ^ Ug du^ ^ Uif, du^ ^ UrU^ ^ UgUfj, cot 9 I dp 

dt ^ dr r d9 r sin 9 d(j) j. ^ ^ gjj^ q 

/ 2 dUr 2 cos 9 dug Uff, 

V sin 9 dS sin^ 9 dd) sin^ i 



1 9 , 2 \ 1 f n \ 1 (^"^4) 

-^ir V ^'■i (smfclug) H ^irf = 0, 



where 



I d f 1 d f . ^d\ 1 52 



A = + , . — sin 



^2 \^ j ^2 gj]^ 9 d9 \ d9 ) sin^ 9 dcj)^ 

is the Laplacian in spherical coordinates and the initial and boundary conditions are not written. 
One can check by direct differentiation that there exists an exact stationary solution to the three- 
dimensional NS equations in spherical coordinates: 

2 

01 ex. 

Ur = G, Ug = — , U^ = 0, p = (3- . „ , (1.3) 

r sm 6* ^ 2r^ sur 9 
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where (a, (3) are arbitrary parameters. The stationary solution (|l.l-{|) describes fluid motion tan- 
gential to a sphere of any given radius r. The stationary flow has two pole singularities at ^ = 
and 6 = TT. The singularities correspond to the source and sink of the velocity vector at the North 
and South poles of the spherical shell 17: the fluid is injected at the North pole from an external 
source and it leaks out at the South pole to an external sink. 

In the limit e ^ 0, the non-stationary three-dimensional fluid flow is confined on a sphere S of 
unit radius parameterized by the polar (latitude) angle 6 and azimuthal (longitude) angle (p, 

S = {{9,(l)), ^ < ^, ^ (/>< 27r} . (1.4) 

Since the velocity vector u and the pressure p in the NS equations Hl.l|) are coupled together by 
the incompressibility constraint V • u = 0, it is difficult to analyze the full set of three-dimensional 
equations. A common approach to simplify the problem is to use the artificial methods such as 
the pressure stabilization and projections [23 . The error estimate of the pressure stabilization 
and projection methods is not however mathematically precise. Instead, we shall use the result of 
the Theorem B in [22], which states that provided the function uo(r, 6',(/>) is smooth enough, the 
strong global solution u(r, 6, cp, t) of the three-dimensional NS equations converges as e ^ to the 
strong unique global solution •v{9,(j),t) of the two-dimensional NS equations on the sphere, where 
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l+e 



v(6',(/),t) = lini - / ru(r, 6',(/),t)dr = (0,fe,f^). 
£ Ji 

The vector v(^, </), t) is interpreted as the average velocity with respect to the radial coordinate r. 
The two-dimensional NS equations on a sphere S in spherical angles [9, (j)) are written explicitly 
as follows |22]: 

dvB dv0 dve 2^0 , f a '"e 2cos9dv^\ 

ot o9 sm 9 ocp ^ d9 \ sm^ 9 sm^ 9 ocp J 

dvA, dvA , VA, dv^ 1 dp f 2cos9 dve va 
+ ^^1^ + — + vev^^oi9 = -—^^ + ly ( Asv^ + —2J1^ - —2, 



dt 89 sin 9 dcf) sini 

ism9ve)+ . ^ J = 0, 



sm U om sm 



sin 9 89 sin ( 

where A5 is the Laplace-Beltrami operator in spherical angles 

1 d f 8\ 1 
sm9— + 



sin 6* 50 y""d9j ' sin2 502- 

Note that no boundary conditions are specified for the vector v(9,(j),t) on sphere S, while the 
initial condition v|j=o = vq on S" is not written. For the purposes of our work, we rewrite the 
two-dimensional NS equations on the sphere S in an equivalent form: 

dvg Vd,uj dq ( Ve 2cos9 dvs\ 

at smfc* 89 \ sixr' 9 sm^ 6* dcf) J 

ot sm 9 sm.9 ocp \ sur 9 ocp sm^ 9 J 

|,.i„,„„) + ^ = 0, (1.7) 
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where g is a static (stagnation) pressure and uj is the vorticity: 

q = P+^{ve+v^) , uj = —{smev^)-—. (1.8) 

The stationary solution H1.3|) corresponds to the exact stationary solution of the two-dimensional 
NS equations (|1.5j) - (|1.7j) on the unit sphere S: 

OL 

ve = V4> = 0, q = P, (1.9) 

where (a, (3) are arbitrary parameters. 

We shall also consider the situation where the external source and sink singularities at ^ = 
and 9 = TT are excluded from the domain of the NS equations p.5|l - (|1.7jl . For instance, we shall 
consider the truncated domain in the form of the spherical layer 

= {(^,0): eo<9 <TT-eo, 0<(P<2tt}, (1.10) 

where < < f • Without loss of generality, the spherical layer Sq is truncated symmetrically at 
the two rings located in the Northern and Southern semi-spheres such that the stationary flow ()1.9|) 
is free of pole singularities in Sq. In other words, without dipping into details on how the fluid flow 
is injected on the sphere and is collected from the sphere in a neighborhood of the North and South 
poles, we will study how the fluid leaks from the Northern semi-sphere to the Southern semi-sphere 
along the spherical layer (|1.1U|) . In this context, the stationary solution H1.9() is interpreted as the 
mass conservation law which is obtained by integrating the free divergence condition (|1.7)) . 

We are interested in spectral stability of the stationary fluid flow H1.9() . In the case of 5 when the 
singularities are included, we prove analytically that the linearized NS equations H1.5() - ()1.7() about 
the stationary solution (|1.9)) are asymptotically stable. In the case of when the singularities 
are excluded, the asymptotical stability of the stationary flow can only be proved for the case 
ly = oo, that is in the limit of zero Reynolds numbers. By using the power series expansions, 
we approximate solutions numerically and show that the stationary flow remains asymptotically 
stable for all Reynolds numbers. 

Our paper is structured as follows. Section 2 introduces the linearization of the two-dimensional 
NS equations (|1.5j) - (|1.7j) at the stationary solution (|1.9j) and discusses boundary conditions for 
the perturbation vector. Analytical results on location of the spectrum of the linearized problem 
are reported in Section 3 for symmetry-breaking ((/>-dependent) perturbations and in Section 4 for 
symmetry-preserving (^-independent) perturbations. Numerical results on computations of eigen- 
values of the linearized problem are described in Section 5 for symmetry-breaking perturbations 
and in Section 6 for symmetry-preserving perturbations. Section 7 discusses applications. 

2 Linearized equations and separation of variables 

Without loss of generality, we consider the stationary solution (|1.9|) with a = 1 and /3 = 0. 
The presence of arbitrary parameters (a, /?) introduces time-independent (neutral) modes of the 
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linearized equations, which we wih also account for in this section. We consider an infinitesimal 
time-dependent perturbations of the stationary flow with a = 1 and /3 = in the form 

ve = -^ + U{e,cl))e^\ v^ = V{e,(l))e^\ q = Q{9,^)e^\ (2.1) 
smt7 

where A G C is a parameter, such that perturbations with Re(A) > imply spectral instability 
of the stationary flow. If Re(A) < for all perturbations, the stationary flow is asymptotically 
stable, while if Re(A) = for some perturbations and Re(A) < for all other perturbations, the 
stationary flow is stable in the sense of Lyapunov. 

By neglecting the quadratic terms of the perturbation, we linearize the NS equations H1.5I) - H1.7|) 
with the expansion to the form: 

.rr dQ U 2cos^5y\ , , 

do \ siW^O snr6' o 



1 (d dU\ 1 dQ IcosOdU V 

— {sme V)- — \+ — = 1^ AsV + — ^TTT - -T- 
sm^ 9 \od oS J sm 6 do \ sm 26 do sur 



|(.i„«t,) + ^=„. (2.4) 

Perturbation terms of the velocity vector must satisfy some boundary conditions in the domains 
5*0 or S. It is naturally to assume that the velocity vector is periodic with respect to the angle (f): 

C/(0,0 + 27r) = f/(0,(/)), V{e,(t) + 2-n) = V{e,<t))- (2.5) 

Therefore, we look for Fourier series solutions of the system H2.2() - H2.4|) : 

u{e,(t^) = Y,Uk{e)e'''\ v{e,^) = Y,yk{e)e'''t Q{e,^) = Y,Qk{ey^'^- (2.6) 

fcez fcGZ fcez 

We also require that the components (f/^, 14) of the velocity vector be square integrable in Sq or 
S with respect to the spherical weight: 



/ {\Uk? + \Vk?)s\nedB <oo, (2.7) 
Jen 



where < Oq < 7r/2. When the domain is the truncated spherical shell 5*0, we require that 
components the velocity vector vanish at the regular end points of the domain: 



Ukieo) = UkiTT- Oo) = VkiOo) = ViTT - Oo) = 0. (2.8) 

The complete sphere S with the singular end points will be considered in the limit ~^ 0. We 
require that the components of the vorticity in (|1.8|1 vanish at the singular end points of the 
domain: 

lim Uk(e) = lim Uk(0) = lim sin 6*14(61) = lim sin 6*14 (6*) = 0. (2.9) 
6^0 e^TT e-^o e^TT 

It will be clear later that separation of variables is different between the cases k = and k ^ 0. 

We say that the correction terms with k = represent symmetry-preserving perturbations of 
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the stationary flow (|2.1|) . while the correction terms with A; 7^ represent symmetry-breaking 
perturbations. 

Case k ^ 0: It follows from the divergence-free condition H2.4j) that one can introduce the stream 
function ^'fc(^) for the velocity vector {Uk,Vk) as follows: 

Uk = ^^k{0), Vk = -%ie). (2.10) 

smfc* 

The system of linearized equations (|2.2jl - (|2.3l) reduces to the coupled ODE system for ^k{G) and 
Qk = ikPk{e): 

^Pk = -^(uAk^k-X'i'k), (2.11) 
ad sm u 

^ -^Pk = (i^Afel-fc - A^-fc) - -^Afc^'fc, (2.12) 



where 



sm U do sm ( 

<f COS0 d k'^ , , 

dO^ sm 6 dO sm^ 

Let <I>fc = Ajt'ffc be a new variable. Then, the variable Pk can be excluded from the system 
()2.11|) - H2.12|) . such that the system reduces to a closed second-order ODE: 

i/Afc$fc = X^k. (2.14) 

Besides the relations (j2.11|) and (|2.12j) between the pressure P^, the stream function and the 
vorticity we note another relation between these components: 

^>fc = Afc^'fc = sin2^AfcPfc. (2.15) 

Due to the boundary conditions (|2.8() and the representation (|2.1U() . the solution ^^(0) for the 
truncated spherical layer Sq is defined on a closed interval 6*0 < 6* < vr — 0o for < < ^/2 
subject to the boundary conditions 

^fe(^o) = ^fc(^o) = ^fc(vr - ^0) = ^'U^ - ^0) = 0. (2.16) 

Since = and 9 = ir are singular points of the system ()2.14|) - H2.15|) when ~^ 0, the solution 
^'fc(0) for the complete sphere S is defined on an open interval < 9 < tt satisfying the boundary 
conditions from (|^ and (I^TTUl) : 

lim ^fc(6l) = lim sin 0^^(61) = lim ^k{0) = lim sm9%{9) = 0. (2.17) 

Case k = 0: It follows from the divergence-free condition H2.4|l that 

a 



U< 







sm y 

where a G R. This solution resembles the neutral eigenmode generated by the arbitrary constant 
a in the stationary solution (|1.9|1 . Since the eigenmode violates the boundary conditions 1)2.16^ on 
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So and has pole singularities on S, we set a = 0. In this case, the first equation (|2.2|) admits a 
solution Qq = f3, where /9 S M. It is also a neutral eigenmode generated by the arbitrary constant 
P in the stationary solution H1.9() . Since it is a trivial eigenmode (the pressure term is defined with 
accuracy to an addition of an arbitrary constant), we can set /3 = 0. 

When a = f3 = 0, the representation for Uq = Qq = matches the previous representation for Uk 
and Qk with A; = 0. Using the representation ()2.1U() . we introduce Vq = —^q{0) and rewrite the 
second equation (|2.3jl as follows: 



^ (i^Ao^-o - A^-o) 



1 



sm ( 



0^0 



0, 



(2.18) 



where Aq is defined by (|2.1,''{j) with k = 0. Letting <I>o = Ao^'o and taking one more derivative in 9, 
one can convert the non-trivial equation (|2.18l) to the previous form H2.14() with k = 0. Therefore, 
all solutions of (|2.18() are also solutions of H2.14() with k = 0, while the converse statement is 
not true. It follows from (|2.8() and 1)2. 9|1 that the stream function ^oi&) satisfies the Neumann 
boundary conditions 

^o(^o) = ^o(^ - 0o) = (2.19) 



in the case of So and the boundary conditions 



limsin6l*f)(6l) = lim sin 6l^''n(6') = 

0^0 e-^TT 



(2.20) 



in the case of S. Stability analysis of the linearized system ()2.14|) - (|2.15() with /c / is developed 
separately from that of the linearized equation (|2.18|) with k = 0. Our main results on eigenvalues 
of the linearized systems H2.14[) - ()2.15() and 1)2. 18(1 are summarized in Table 1. The remainder of 
this article is devoted to the proofs and numerical verifications of results described in Table 1. 



Index k 


Viscosity i' 


Cut-off 00 


eigenvalues 


results 


k^O 


< u < oo 


00 = 


real negative 


Proposition 2 


k^O 


V = CO 


< 00 < f 


real negative 


Propositions 3 and 4 




< < oo 


< 00 < f 


real or complex 


Section 5 


A; = 


Q <u < oo 


00 = 


real negative or absent 


Proposition 8 


k = 


Q <u < oo 


< 00 < f 


real negative 


Propositions 9 and 10 


k = 


< < oo 


< 00 < f 


real negative 


Section 6 



Table 1: Summary of main results. 



3 Stability analysis for k ^ 

We rewrite the coupled system ()2.14() - (|2.15() for {^k, ^k) by using the variable x = cos 9: 

ifc*fc = ^fc, Lk'^k + e'^'k = f^^k, (3.1) 
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where e = 1/v is the Reynolds number of the basic flow, ^ = \/v is a rescaled eigenvalue, and Lk 
is the Sturm-Liouville operator for associated Legendre functions 



Lj^ — — 
dx 



ax 



1 — X 



2 • 



(3.2) 



The system (|3.1|1 is defined on the symmetric interval —xq < x < xq, where xo = cos^o- The 
spherical layer 5*0 corresponds to the case < xq < 1, while the complete sphere S corresponds to 
the limit xq — > 1. In the latter case, the interval x £ [—1,1] connects two singular points x = ±1 
of the Sturm-Liouville operator (|3.2() . The case e = corresponds to the infinitely viscous fiuid, 
while the case e = oo corresponds to the inviscous fluid. 

Using the representation 1)2. 1U() and the transformation x = cos e with %{e) = -Vl-x^%{x), 
we rewrite the condition 1)2. 7() as the norm on function space TCk, which is used throughout our 
work: 



1^- 



,|2 



XQ 



XQ 



(1 - x')\%{x)\' + ^\^,{^ 

1 — 



dx < oo. (3.3) 



We shall denote Tikii—XQ, xq]) when < xq < 1 and 7ik{[—^, 1]) when xq = 1. When < xq < 1, 
the linearized system (|3.1|) is defined on function space 

Xo = {^-fc G Wfc([-xo,xo]) : ^k{±xo) = ^'k{±xo)=0} , (3.4) 

where the boundary conditions ()2.16|) are taken into account. When xq = 1, the linearized system 
()3.1|) is defined in function space 

X = l^-fc G Hkii-l, 1]) : lim ^'^(x) = liin (1 - x^)%{x) = o] , (3.5) 

x^±l x-^±l J 

where the boundary conditions ()2.17|) are taken into account. We note that the boundary conditions 
in the definition of X are redundant, since the norm ()3.3|) is finite on x G [—1,1] only if the boundary 
conditions in ()3.5() are satisfied. Nevertheless, we write these redundant boundary conditions 
according to the standard formalism of the singular Sturm-Lioville problems j2Uj . 

The Sturm-Liouville operator in ()3.2|) is self-adjoint with respect to the boundary conditions in 
Xq and X, such that (^'fc,Lfc^'fe) = — ||^fc||>^^ < is finite and real-valued for ^'^ G Tlk{[—xo, xq]). 
Therefore, the kernel of is empty in Xq and X. Because the smallest eigenvalue of L/. is bounded 
away zero, the operator is invertible and range(Lfc) is dense in the space of square integrable 
functions on x G [— Xo,xo] for any < xq < 1. Therefore, as it follows from the first equation of 
the system (j3.1|) . the component $ G range(Lfc) is square integrable on x G [— xo,xo] but it does 
not satisfy any specific boundary conditions at the end points x = ixg. 

The eigenvalue problem 1)3. 1(1 in Xq and X has two continuous parameters < xq < 1 and e > and 
one integer parameter k G Z\{0}, while is the eigenvalue-eigenfunction pair that defines 

spectral stability of the stationary fiow. The following results characterize the spectrum of the 
eigenvalue problem in the cases: (i) xq = 1 and e > 0; (ii) < xq < 1 and e = 0; and (iii) in the 
limit Xq — > 1 when e = 0. Based on these results, we prove the following theorem: 
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Theorem 1 When xo = 1 and e>0 or 0<xo<l and e = 0, the stationary flow is 
asymptotically stable with respect to symmetry-breaking perturbations in the sense that the spectrum 
of the linearized problem in Xq or X consists of a set of isolated eigenvalues /i of finite 

multiplicities, where ^ G M_ is bounded away from zero. 

The proof of theorem consists of the proofs of three individual propositions. 

Proposition 2 A complete spectrum of the eigenvalue problem h'j. 1)) with xq = 1 and e > in X 
consists of simple isolated eigenvalues at jj, = Hn, 

= -Sn(Sn + 1), Sn = CF + n, (3.6) 

where a = y^fc^ + > and n> is integer. 
Proof. Let ^ = —s{s + I) and 

^>^(^)={\^y'^(^)- (3-7) 

The second equation of the system 1)3. 1|) transforms to the associated Legendre equation 

-^ip + sis + l)ip = 0, -l<x<l, (3.8) 

where a = \/k'^ + > 0. Since the hnear ODE (|3.8j) has no singular points on — 1 < x < 1, 
there exists a set of two linearly independent twice continuously differentiable solutions in any 
compact subset ofxS(— 1,1) Singularity analysis of the ODE (|3.8() as x — > ±1 shows that the 
solution (p{x) either have a singular (unbounded) behavior like (1 =F x)~~"^'^ as x ^ ±1 or a regular 
(vanishing) behavior like (1 =F x)'^' ±1. 

Let <p{x) be a regular solution of (|3.8|) on x G [—1, 1], such that tf{x) ~ (1 =F x^^^ and ^k{x) ~ 
(l^x)^''^^'^'^^'^ as X ^ ±1. Since the Sturm-Liouville operator is invertible on <I>fc G L'^([—l, 1]) 
for k ^ 0, the first equation of the system 1)3. Ijl admits a solution ^a:(x) that behaves like (1 =F 
ri.'ji±e/^+cT/2 2; _> Since ±e + Ve^ + Ak"^ > for any A; G Z and e > 0, the function ^k{x) is 
bounded on x G [—1,1] such that G ^^([—1)1]) while the function ^'a,.(x) belongs to the function 
space X in (|3.5j) . Therefore, if (p{x) is a regular solution of (|3.8|1 . then ^'fc(x) is an eigenfunction 
of the eigenvalue problem (|3.1() in X. 

Let (p{x) be a singular solution of (|3.8|) . such that ip{x) ~ (1 =p x)^'^/^ and <I>fc(x) ~ (1 =F x)^^/^~'^/^ 
in at least one limit x — > ±1. Since ite — \/e^ + < —2\k\ < —2 for k ^ and e > 0, the 
function does not belong to L^([— 1, 1]) and hence ^'fc(x) can not be in X. By Theorem 10 on 
p. 1441 in jlUj . the essential spectrum of the formally self-adjoint operator (|3.8() is void. Therefore, 
the complete spectrum of the linearized system (|3.1|) in X consists of isolated eigenvalues /i, which 
correspond to regular solutions ^p{x) of the associated Legendre equation 1)3. 8() . 



d_ 

dx 



^dif 
dx 
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Let ip{x) be a regular solution of ()3.8() and write (p{x) = (1 — where F{x) is bounded 

as X — > ±1. This substitution transforms the associated Legendre equation (|:-{.8j) to the hypergeo- 
metric equation 

z{l - z)F"{z) + (7 - (a + /3 + l)z) F'{z) - a/3F{z) = 0, (3.9) 

where 

1 — X 

z = — - — , a = a — s, (3 = a + s + l, 7 = 0" + !. (3.10) 

The only solution of the ODE H3.9|) which is bounded as j; — > 1 (z ^ 0) is the hypergeometric 
function F(z;a,/3, 7), which admits the power series at 2; = (see 9.100 on p. 995 in [T^l: 

Fiz;a,P,^) = l + —z + ^-^-—^z + ... (3.11) 

The hypergeometric series (|3.11|) converges for \z\ < 1 but it diverges as z ^ 1 since a + /9 — 7 = 
£7 > unless the truncation of the power series to a polynomial in z occurs (see 9.101-9.102 on p. 
995 in The latter case is the only case when the solution of the ODE 1)3. 9(1 is bounded in both 

limits X ^ 1 (z — > 0) and x — > — 1 (z — > 1). It is easy to see that the truncation occurs when either 
a = —n or P = —m with non-negative integers n and m. The two cases are in fact equivalent to 
each other since fi = —s{s+l) = {a—a){P—a) and a+P = l+2a. Let a = —n, such that s = a+n, 
/? = 2(T + n + 1 and 7 = o" + 1. In this case, the function F (z; — n, n + 1 + 2a, 1 + o") = Fn{x) is a 
polynomial of degree n, e.g. 

' 2(1 + 0-) ' 2(1 + 0-) ' ^ ' 

while the simple eigenvalues /i = /u„ are given by the expression 1)3. 6|) . When o" = 0, polynomials 
Fn{x) coincide with the Legendre polynomials Pn{x) in 8.91 on p. 973 of ^3]. ■ 

Proposition 3 A complete spectrum of the eigenvalue problem h'j.l]) with < xq < 1 and e = in 
Xq consists of isolated eigenvalues fi, which are (i) real and strictly negative and (ii) either simple 
or double with linearly independent eigenf unctions. 

Proof. We first show that no zero eigenvalue = exists in the eigenvalue problem ()3.1() with 
< xo < 1 and e = in Xq. Let ^'^(x) be a C^([— xo,xo]) solution of the fourth-order ODE 
L\^k = in function space Xq. Then, 

{^k,Ll^k) = (1 - x^) [^kiLk'^k)' - ^U^fc^fc)] \l=-\ + {Lk^k^Lk^k) = {Lk^k^Lk^k)- 

Therefore, ^'fe(x) is in fact the solution of the second-order ODE Lk^k = 0. The boundary 
conditions in Xq admit the only solution "^k^x) = 0, such that the eigenvalue problem (|3.1|) 
contains no eigenvalue /i = in Xq. 

When /i 7^ and e = 0, the system (|3.1() admits a general solution in the form 

^k{x) = ^ + Hx), ^k{x)=<P{x), 
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where ^l^{x) and (j){x) are general solutions of the homogeneous second-order ODEs 

Lktp = 0, Lk<p = ficj). 

Since the operator is invariant with respect to the inversion symmetry x i— > — x, each homo- 
geneous second-order ODE has linearly independent symmetric (even) and anti-symmetric (odd) 
solutions denoted by subscripts + and — respectively. Therefore, we obtain the decomposition 

^k{x) = d+^i^ + c+V^+(x) + d_^^+c_V-(x), 

^kix) = (i+(/)+(x) + ci_0_(x), 

where (c+, c_, d+, (i_) are constants and the functions (j)±{x) and 'il}±{x) are uniquely normalized 
by the initial values at x = (e.g. 0+(O) = 1, 0+(O) = and 0_(O) = 0, (j)'_{0) = 1). We note that 
either 'ip±{xQ) 7^ or ^^(xq) / (since ip±{x) = otherwise). By using the boundary conditions 
in (|3.4|1 . we decompose the boundary- value problems into two uncoupled systems with 

d±4>±{xo) + fJ'C±ip±{xo) = 0, d±cl)'^{xo) + /UC±V±(a;o) = 0, 

such that a non-zero solution for (c+, c_, (i_) exists provided 

(l>±{xo)'>p±{xo) = (l)±{xo)Tp'^{xo). 

The functions ij±{x) are independent of /i, while (j)±{x) depend on /u. We have thus obtained that 
the functions <j)±{x) solve the cZosefi eigenvalue problem 

Lfc(/)± = fi4>±, -xo <x<xq, (3.13) 

defined on the function space 

Ho = {(j)± enk{[-xo,xo\) : ^p±{xo)(p'±{xo) - i^'^{xo)(/)±{xo) = 0, (/)±(-x) = ±0±(x)} . (3.14) 

The //-independent boundary values in ()3.14|) are Robin boundary conditions when ^'^{xo) and 
iP'^(xq) are both non-zero, Dirichlet boundary conditions when tp±{xo) = and Neumann boundary 
conditions when ^p'j.{xo) = 0. The associated Legendre operator Lk is self-adjoint in Hq with respect 
to any of these boundary conditions [20]. Therefore, all eigenvalues n of the eigenvalue problem 
()3.13|) in Hq are real- valued and isolated, while the corresponding eigenfunctions 4'±{x) are real- 
valued. Moreover, all eigenvalues of 1)3. 13(1 are simple since the Wronskian of any two solutions of 
p.l3|l with boundary conditions in ()3.14|) is zero. Since ^'fc G Xq and ^k £ Hq, we obtain that 

((/., (p) = {Lk^k, </>) = (^fc, Lk(l)) = /i(^'fc, 0) = ((/>, 0) + iMiil), 0), 
such that =0 for fJ- ^ 0. By using the above identity, we obtain that 

^(0, cj)) = i^k, 0) = (^fe, Lk^k) = -W'^kWn, < 0, (3.15) 

such that /u < for each eigenvalue with "ifk ^ and (/> / 0. By construction, eigenvalues are 
at most double. The case of double eigenvalues corresponds to the situation when the eigenvalue 
problems H3.13() - (|3.14p admit two linearly independent (even and odd) eigenfunctions for the same 
value of ■ 
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Proposition 4 Let {^n}n>Q be isolated eigenvalues of the eigenvalue problem i,V. 1\) in Xq with 
< xo < 1 and e = ordered as 

/"o > /^i > ••• > /in > ••• 

Then, 

lim /x„ = -Sn(s„ + 1), Sn = \k\+n, 

Xo^l 

where n > 0. 

Proof. Consider even and odd solutions of the second-order ODE Lkip± = in the hmit xq ^ 1. 
Since the kernel of admits no eigenfunctions in Ti^ for A; 7^ and < xq < 1, the solutions 
V'±(xo) must diverge as xq 1. Singularity analysis as x — > ±1 suggests that the solution '4)±{x) 
grows like (1 =F x)"''^'/^ as x ^ ±1, such that lim V'±(2;o)/V'±(^o) = 0. Therefore, eigenfunctions 

a:o— >1 

(t>±{x) of the auxiliary eigenvalue problem H3.13(l for < xq < 1 satisfy in the limit xq ^ 1 the 
singular eigenvalue problem 

Lk4>± = l^4'±i — 1 < X < 1 (3.16) 

defined on the function space 

H={<P±enk{[-l,l]): lim </.±(x) = hm (l-xV±(^) =o|. (3.17) 

Again, the boundary conditions in H are redundant due to convergence of the integral in 7ik{[—l, 1]). 
A complete spectrum of the eigenvalue problem ()3.16() ~ ()3.17() is constructed in the proof of Proposi- 
tion[2| eigenvalues are given by (|3.6|) with e = and eigenfunctions are (f>±{x) = (1 — x^)l^l/^-F„(x), 
where Fn{x) are associated Legendre polynomials ()3.12() with a = \k\. Convergence and unique- 
ness of continuations from eigenvalues of H3.13() in Hq for xq < 1 to eigenvalues of 1)3. 16(1 in H 
for Xq = 1 is proved in two steps. Theorem 5.3 of j24j guarantees convergence and uniqueness 
of continuations from the singular Sturm-Liouville problem (j3.16|) in H to the regular Dirichlet 
problem for the Sturm-Liouville operator H3.13() on — xq < x < xq. The Dirichlet problem is gener- 
ally different from the Robin boundary- value problem in Hq by the terms 'tp±{xQ)(j)'^{zizXo) /^l^'j.{xQ) 
in the boundary conditions in Hq. However, these terms are small in the limit xq — > 1. Unique 
continuation of simple eigenvalues of the Dirichlet problem to the simple eigenvalues of the Robin 
problem (separately for 0+(x) and 0_(x)) follows by standard perturbation theory of eigenvalues 
of self-adjoint Sturm-Liouville operators in Lemma VIII 1.24 of [15) . ■ 

Remark 5 Theorem ^ does not cover the case < xq < 1 and e > 0. Eigenvalues of the linearized 
problem (13.1)1 in this case will be computed in Section 5 numerically. 

4 Stability analysis for k = 

We rewrite the linearized equation ()2.18)) in the variable x = cos 9: 

Lq^q = $0, ^0 + TT^^o = l^%- (4.1) 
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where Lq is the Sturm-Liouville operator for Legendre functions 

d 



Lo = - 

dx 



(4.2) 



and ^q{x) is introduced similarly to the system Incorporating the boundary conditions 

()2.19|) and H2.2U() in new variables, we introduce the function spaces Xq and X for the eigenvalue 
problem (|4.1|1 . When < xq < 1, the function space Xq is 

Xq = {^0 G 7^o([-^o, ^o]) : ^'o(±^o) = 0} . (4.3) 
When xq = \, the function space X is 



X = 1^-0 G Ho([-l, 1]) : ^nm^(l - x')^'^{x) = o| 



(4.4) 



where the boundary conditions are redundant due to convergence of the integral in 7iQ{[—\^ 1]). No 
boundary conditions on ^'0(2;) are set at x = ixo- Moreover, the system (|4.1|) defines the function 
^q{x) up to an arbitrary additive constant. Therefore, the constant function ^q{x) = const is 
always an eigenfunction of the system 1)4. 1() with ^q{x) = 0. 

Lemma 6 The eigenvalue fi = of the linearized system i4.1\ ) in either Xq or X is algebraically 
and geometrically simple. 

Proof. Integrating the first equation in the system (|4.1j) on x G [— xo,xo] for "i^oix) in either Xq 
or X, we obtain the Fredholm Alternative condition 

^o(x)dx = 0, (4.5) 

Xo 

where < xq < 1- Integrating the second equation in the system ()4.1|) . we obtain a general solution 
for n = 0: 

/l-x^'/2 

= (tt^^ 

where cq is constant. Since ^o{x) does not satisfy the Fredholm Alternative condition ()4.5() . we 
have to set cq = 0. Then, ^'o(x) satisfies the second-order ODE Lq^o = 0, which admits only one 
eigenfunction ^'o(x) = const in either Xq or X. Similarly one can prove that the Jordan block of 
the zero eigenvalue with the eigenfunction ^'o(x) = const and ^o{x) = is of the length one. ■ 

We will extend results of Section 3 to the linearized problem 1)4.11) with 7^ in Xq and X. 
Neglecting the only zero eigenvalue /x = with the trivial eigenfunction ^'o(x) = const, we prove 
the following theorem. 

Theorem 7 The stationary flow lil.9\) is asymptotically stable with respect to symmetry-preserving 
perturbations in the sense that all eigenvalues /i (excluding the trivial zero) of the linearized problem 
i4.1[ ) with < Xo < 1 and e >0 in Xq or X are real and strictly negative. 
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In order to develop analysis of eigenvalues for ^ 7^ 0, we shall use two equivalent reformulations of 
the third-order ODE system ()4.H) as the second-order eigenvalue problems associated with formally 
self-adjoint operators. In the first reformulation, we exclude /i\I'g(a;) from the system ()4.1|) and 
find a closed equation for $0(2;), 

Lo<^o + e^o = Ai^o- (4.6) 



By introducing new dependent variable ip{x) via 



e/i 



(4.7) 



1 + x 

the linearized equation (|4.fcij) is transformed to the self-adjoint form given by the associated Leg- 
endre equation 

rl ^ r,^ dipl , , , 

-xo < X < xq. (4.8, 



d_ 

dx 



{l-x')f- 
dx 



-if = fnp, 



4{l-x\ 

By using the second equation of the system H4.1|) . we obtain the first-order ODE for the function 
^o{x): 



H%ix) 



1 



x 



1 + x 



e/4 



dip e 
d^ ^ 2(1 - 



x2) 



(4.9) 



While the linearized equation (|4.6() coincides with the second equation of the system H3.1|) for 
k = 0, the present role of this equation is different. In order to find ^o{x) from a solution ^o{x) 
of the closed equation 1)4. 6|) . we can solve the first-order ODE ()4.9() in either Xq or X with /U 7^ 0. 
Therefore, as opposed to the case A; 7^ 0, we do not have to solve the first equation of the system 
()4.1|) and the Fredholm Alternative condition ()4.5() can be ignored in this approach. 

In the second reformulation of the third-order ODE system (|4.1I) . we introduce a new dependent 
variable xi^) via 



1 



x 



1 + x 



e/4 



(4.10) 



By using the first equation of the system (|4.H) . we express the function ^q{x) in terms of xi^): 



$n(x) 



1 



./4 



:dx 



e + 2x 



, 1 + X y V ' " dx 2\f\- X- 
The second equation of the system (|4.1|) transforms then to the self-adjoint form: 



(4.11) 



d_ 

dx 



dx 



+ 4 + 4ex 
4(1 - x2 



-X = f^X, 



-Xo < X < Xq. 



(4.12) 



Although the second-order ODE 1)4. 12() is more complicated than the associated Legendre equation 
()4.8|) . the eigenfunction x{x) is related to the function ^'o(x) better than the eigenfunction ip{x). 
In particular, when xq = 1 and ^'o € X, the eigenfunction x(x) satisfies the conditions: 



1 



X 



1 + x 



e/2 



X^ix)dx < 00, 



lim 



±1 V 1 + X 



e/4 



■\/l — x'^xix) = 0. 



(4.13) 
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When < xo < 1 and S Xq, the eigenfunction x{^) is any classical solution of the second-order 

ODE (|4.12() on X £ [—xo,xo] with the Dirichlet boundary conditions xi^^o) = 0- There exists a 

pair of Darboux-Backlund transformations between the Sturm-Liouville problems (|4.8() and H4.12() : 

. £ _|_ 2x 

^p{x) = VI - x^x'{x) - ^ ^- (4.14) 



2VI-X 



fix{x) = \/l - x^y'(x) + ^ ^__ =y(x), (4.15) 



where 7^ is assumed. By the Friedrichs' theorems (see, e.g. Theorem 10 on p. 1441 or Theorem 
67 on p. 1501 of ^0]), the essential spectrum of the formally self-adjoint operators 1)4.8(1 and (|4.12|) 
is void. Therefore, the spectrum of these operators consists of a sequence of isolated eigenvalues 
of finite multiplicities, which we identify in three individual propositions. 

Proposition 8 A complete spectrum of the eigenvalue problem J^. j| ) with xq = 1 and < e < 2 

in X consists of simple isolated eigenvalues at = fin, where 

/^„ = -n(n + l), n>0. (4.16) 
No non-zero eigenvalues of the eigenvalue problem i4.1\ ) with xq = 1 and e > 2 exists in X. 

Proof. Let = —s{s + 1) 7^ and ^{x) = (1 — x'^Y^'^F{x) and consider the associated Legendre 
equation (|4.8() with xq = 1. Then, the function F{x) satisfies the hypergeometric equation (|3.9|) 
under parametrization (|3.10|) with a = e/2. In order to identify solutions F{x) of the hyperge- 
ometric equations in the function space ^'q G X, we shall rewrite the relation (|4.9|) as follows: 

-s{s + l)^'^{x) = (1 - xy/^ (f'(x) + ^^^^F(x)) . (4.17) 

Also recall that $o(x) = (1 - x)^/^F(x). When e = 0, we find that ^o{x) = F{x) and ^q{x) = 
— j(j^F[x) + const, such that £ X \i and only if F{x) G X. The only set of eigenfunctions 
of the Legendre equation (|4.8() with e = in X is the set of Legendre polynomials F = Pn{x) for 
s = n with n > (see 8.91 on p. 973 in 13 ). This set corresponds to the eigenvalues (|4.16|) . 
Although the zero eigenvalue (s = n = 0) is excluded from the approach above, it is still added to 
the spectrum by Lemma IHl 

When e > 0, the eigenfunction ^'o(x) belongs to X only if F{x) has a regular behavior as x = 1 
[z = 0). The only solution of the hypergeometric equation (|3.9() which is bounded as x ^ 1 is 
the hypergeometric function F(z;a,/3, 7). (Indeed, by 9.153 on p. 1001 of [T3|, the other linearly 
independent solution F{x) has a singular behavior like F{x) ~ (1 — x)"*^/^ as x — > 1, which results 
in the divergence ^o(x) ~ (1 — x)^"*^ as x ^ 1, such that ^'o ^ X.) By the identity 9.131 on 
p. 998 of US]) the hypergeometric function F(2;;a,/3, 7) admits the following behavior at the other 
singular point x = — 1 (2: = 1): 

Th)Th -a- 0) 

F{z- a, 13, 7) = ^ ^^(1 - z; a, /3, a + /? - 7 + 1) 

r(7 - a)r(7 - /?) 

+(1 _ ^)7-a-/^ r(7)r(a + /3 -7) ^(^ _ ^. ^ _ ^ _ ^ _ „ _ ^ + 1), (4.18) 
T{a)T{l3) 
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where T{z) is the Gamma function and 

1 — X e ^ e ^ e ^ 

Since a + /? — 7 + 1= 7 and 7 — a — /?+l = l — |, the relation ()4.18p can be used only for e < 2 
(the hypergeometric function F(z;a,/3, 7) diverges for 7 = —n with n > integer). 

It follows from H4.17() that the first term in H4.18() leads the singular behavior of ^'o(x) ~ (1 + x)~^ 
as X — > —1 (z ^ 1) if e 7^ 0, while the second term in H4.18() leads to the singular behavior 
^'g(x) ~ (1 + x)~'^/2 as X — > —1 (z — > 1) if /i 7^ 0. Therefore, the eigenfunction ^'o(x) belongs 
to X only if the first term in (|4.18j) is removed which is only possible if7 — a = l + s = —n or 
7 — /3 = —s = —m with integers n, m > 0. Both choices define the same set of eigenvalues ()4.16() 
in the parametrization /i = —s{s + 1). Using another identity 9.131 on p. 998 of jl3j . 

F{z; a, P, 7) = (1 - zr'-"~^F{z; 7 - «, 7 - /3, 7), (4.19) 

we set s = —1 — n with n > 1, such that 

F (z; I + 1 + n, I - n, I + 1) = (1 - z)-'/^F (z; -n, n + 1, | + l) , 

where F (z; — n, n + 1, 1 + e/2) = Fn{x) is a polynomial of degree n, e.g. 

- _ _x + a ~ _ 3x2 + 3ctx + 0-2 - 1 i^^3 _^ 25o-x2 + (Gcj^ - 9)x + a{a'^ - 4) 

°" ' '"TT^' (l + a)(2 + a) ' (l + a)(2 + c7)(3 + cT) ' 

with cj = e/2. When e = (cr = 0), polynomials F„ coincide with Legendre polynomials Pn{x) in 
8.91 on p. 973 of jl3j . The zero eigenvalue {n = 0) is excluded from the construction but added 
to the spectrum by LemmalHI When e < 2, the resulting eigenfunction \l!'o(x) belongs to X. Since 
^'q(x) ~ (1 + x)~'^/2 as X — > —1, the resulting eigenfunction \I'o(x) does not belong to X for e > 2. 

We shall prove that no non-zero eigenvalues exist in X for e > 2. Using the identity (|4.19|) . we 
transform the solution F{x) to the equivalent form F{x) = (1 + x)~^^'^F{x), where F{x) satisfies 
the hypergeometric equation 1)3. 9|) with new parameters 

1 — x_ ^ ~ ^ _ e 

2; = —^, Q = 7-a = l + s, /? = 7-/3 = -s, 7 = 7 = - + 1. 

Up to a constant factor, F{x) is represented by the hypergeometric function F{z] 1 + s, — s, l + e/2). 
It follows from the ODE 1)4. 17() that the eigenfunction *I'o(x) is related to F{x) by 

-s{s + l)^',{x) = F'{x). 

Since a + /3 — 7 = — e/2 < for e > 0, the hypergeometric series for the function F{z; 1 + s, — s, 1 + 
e/2) converges absolutely on the entire interval x G [—1, 1] {z S [0, 1]) (see 9.102 on p. 995 of ^13j). 
Therefore, F{x) G on x G [-1, 1] and F'{-1) is well-defined. We shall prove that F'{-1) / 
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for any s 7^ and e > 2. It follows from the hyper geometric equation H3.9() with (q,/3,7) at z = 1 
that 

.(. + l)F(-l) + (|-l) F'(-1)=0. 

If F'{—1) = 0, then F{—1) = for any s 7^ and e > 2, and the only regular solution of the 
hypergeometric equation H3.9() is F{x) = 0. Therefore, F'(—l) ^ 0, and therefore, ^0 ^ X for 
e > 2. ■ 

Proposition 9 A complete spectrum of the eigenvalue problem J^. wii/i < xq < 1 and e > 

in Xq consists of simple isolated eigenvalues ^ with /i G M_ . 

Proof. When ^'o G Xq, the eigenfunction (/^(a;) of the associated Legendre equation (|4.8|) satisfies 
the Robin boundary conditions 

2{l-xl)ip'{±xo) + eip{±xo) = 0, 

while the eigenfunction xi^) of the second-order ODE (|4.12j) satisfies the Dirichlet boundary 
conditions xi^^o) = 0. Each eigenvalue problem is self-adjoint with respect to these boundary 
conditions [20]. Therefore, all eigenvalues /i of the regular boundary- value problems are real- valued 
and isolated. Moreover, these eigenvalues are negative due to the Green's identity pUj : 

^ ip\x)dx = - {l-x^){ip'{x)f dx-- -^^dx<0. (4.20) 

These eigenvalues are also simple, since the Wronskian of any two solutions with the Robin or 
Dirichlet boundary conditions is zero. ■ 

Proposition 10 Let {^n}n>o be isolated simple eigenvalues of the eigenvalue problem i4-l^ with 
< xq < 1 and Dirichlet boundary conditions x{^xq) = ordered as 

> /io > Ail > ••• > Ain > ••• 

Then, lim fin = —Sn{sn + 1), where 

s„ = 1 + n, /or < e < 2 

with n > 0. 

Proof. Singularity analysis of the second-order ODE (|4.12|1 shows that the solution x{x) behaves 
as 

X ^ c+(l - + c+(l - as X ^ 1 

and 

X ^ q (1 + + C2 (1 + as X ^ -1 



and 



e 

— ^ n 
2 



-|- n, for e > 2 
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The ODE (|4.12|) admits a bounded (regular) solution x{^) on x £ [— only if the singular 
components are removed. This leads to the constraints = and either = for < e < 2 or 
C2 = for e > 2. It is explained in Proposition |H1 that the set C2 = and cj~ = for < e < 2 
is equivalent to s = m with m > 0, when the first term in the relation ()4.18|) is removed and 
the hypergeometric function F(z;a,/3, 7) is a polynomial. Note that the zero eigenvalue s = 
(m = 0) of the problem 1)4.8(1 is excluded from the spectrum of the problem ()4.12|1 . such that 
s = Sn = ^ -\- n with n > 0. On the other hand, the set = and = for e > 2 is equivalent 
to s = Sn = e/2 + n with n > 0, when the second term in the relation (|4.18() is removed and 
the hypergeometric function F{z; a, (3, 7) is a polynomial. Note that the first Darboux-Backlund 
transformation ()4.14|) implies that if = 0, then 



^^c+(l 



xY^^, as X 



1 



and 



c5"(l + xy/^ + c^(l + x)-'/^ as X ^ -1. 



Recah that ip{x) = (1 - x^)^/''F(x). When q 



[1 + x) 



-e/2 



3jS X 



(0 < e < 2), F{x) is singular like F{x) 

e > 2), F{x) 



-1 in accordance with the relation ()4.19() . When = 



is bounded as x — > —1. The marginal case e = 2 corresponds to the case when xi^) has a 
bounded and logarithmically growing components as x — > —1. The logarithmic growth is excluded 
if s„ = 1 + n with n > 1, which is the border between the two spectra at e = 2. 

When s = Sn and e 7^ 2, the eigenfunction xi^) of the formally self-adjoint problem (|4.12() satisfies 
the Dirichlet boundary conditions lim^;^-!-! ^(2;) = 0. When e = 2, the eigenfunction xi^) is 
bounded at x = — 1 and zero at x = 1. In either case, convergence and uniqueness of continuations 
from eigenvalues of the regular Dirichlet problem ()4.12|) with xq < 1 to eigenvalues of the singular 
boundary- value problem (|4.12|) with xq = 1 is proved by Theorem 5.3 of j24j . ■ 

Remark 11 Bounded (for e = 2) and decaying (for e > 2) eigenfunctions xi^) of the self-adjoint 
problem (|4.12() with xq = 1 for eigenvalues n = — s„(s„-|-l) with s„ = e/2+n violate the conditions 
()4.13|) . Indeed, one can check that the limit in ()4.13|) as x — > — 1 is non-zero (proportional to c^) 

hence diverges. Therefore, the eigenvalues of the self-adjoint problem 



e/4 



and the integral in ([4.1.'- 
()4.12|) for e > 2 do not correspond to eigenvalues of the original problem (|4.1() in space £ X , in 
agreement with Proposition |H1 We also note that if one consider a generalized conditions for xi^) 
with 

1 — X 



lim 



1 + X 



\/l — x^x{^) 



< 00, 



(4.21) 



the spectrum of the self-adjoint problem (|4.12() is not defined since C2 ^ and the eigenvalue 
problem is not complete. 



Remark 12 Theorem |3 covers the entire parameter domain < xq < 1 and e > 0. However, 
there is an interesting problem with convergence of eigenvalues of the associated Legendre equation 
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(|4.8|) in the limit xq — > 1. While the eigenvalues with < xq < 1 are expected to converge to the 
eigenvalues in H4.16|) for < e < 2, no eigenvalues with the eigenfunctions ^ G X exist for e > 2. 
Convergence of eigenvalues of the linearized problem ()4.1() as xq — > 1 will be computed in Section 
6 numerically. 



5 Numerical computations of eigenvalues for A; 7^ 

In order to illustrate distribution of eigenvalues in Propositions [2 |31 and \^ and to investigate 
eigenvalues in the domain < xq < 1 and e > in Remark [21 we develop a numerical method 
based on power series expansions. Since x = is an ordinary point and x = ±1 are regular singular 
points of the system the power series expansions of the functions ^'fc(x) and ^k{x) in powers 

of X converge uniformly and absolutely for |x| < 1. The numerical method is based on truncation 
of the power series. 

Let ^ G C be parameterized by /i = —s{s + 1), s G C. Due to the symmetry, it is sufficient to 
consider the domain {s G C : Re(s) > — |}. The stability domain Re(/x) < corresponds to the 
domain 

|s G C : |Im(s)| < VRe(s)(Re(s) + 1), Re(s) > o| . (5.1) 
Consider the power series with separated even and odd terms: 



^kix) = ^c^x^^+^d^x^^+i, (5.2) 

m>0 m>0 

<^k{x) = ^a^x2-+^5™x2-+i, (5.3) 

where the starting coefficients (ao, 601 cq, do) are parameters. Substituting ()5.3|) into the second 
equation of the system (|3.1j) we find that (ai, bi) are defined separately as 

_ (fc^ - s(l + s))ao - ebp 

+ 2 - s{s + l))bo - e2ai 
h = , (5.5) 

while the coefficients {am,bm}m>2 are defined uniquely from the recurrence equations: 
(^2 - s{s + 1) + 2(2m + 2f)am+i + {s{s + 1) - 2m(2m + l))o™ 



(2m + 4)(2m + 3) 
-e(2m + 3)6m+i + e(2m + 1)6^ 
(2m + 4)(2m + 3) ' 
£^ - s(s + 1) + 2(2m + 3)^)6^+1 + (s(s + 1) - (2m + 2)(2m + l))b„ 

(2m + 5)(2m + 4) 
-e{2m + 4)am+2 + e{2m + 2)0^+1 



(5.6) 



(5.7) 



(2m + 5)(2m + 4) 

We note that the initial equations H5.4I) - H5.5|) follow from the recurrence equations ()5.6|) - H5.7() for 
m = — 1 with a_i = 6_i = 0. 
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Substituting (|5.2|) into the first equation of the system we find that the coefficients {cm, dm}m>2 
are defined from the coefficients {am, &m}m>o by the recurrence equations: 

_ (fc^ + 2(2m + 2)2)cm+i - 2m(2m + l)cm + flm+i - 

Cm+2 — 77^ , .N/r, —:r\ ^^-^ 

(2m + 4j (2m + 3j 
(/c2 + 2(2m + 3)2)(i„+i-(2m + 2)(2m + l)(i™ + 6™+i-6^ 

'^'"+2 = 7^ — rru^ — rT{ • ^-9) 

{2m + 5) (2m + 4) 

The initial equations for (ci, di) follow from the recurrence equations H5.8() - ()5.9|) for m = — 1 with 
a_i = 6_i = c_i = = 0. 

The boundary conditions in H3.4I) lead to the equations 

J]c„x2- = 0, ^d„x2™ = 0, j;(2m)c^x^™ = 0, (2m + = 0. (5.10) 

m>0 m>0 m>0 m>0 

There exists a linear map from (oq, cq, do) £ parametrized by s G C to the sequence 
{am,bm,Cm,dm}meN- Therefore, the boundary conditions ()5.10() are equivalent to the homoge- 
neous system Afc(s)x = 0, where x = (oq, bo, cq, do)'^ S and Ak{s) is a 4-by-4 matrix computed 
from the entries of (|5.10|) . The matrix ^^(s) depends on s G C and A: G N, as well as parameters 
xq and e. If the power series are truncated at the M-th term, the matrix Ak{s) depends also on M. 
Eigenvalues fi = —s{s + 1) of the system in (|3.4|) are equivalent to roots s of the determinant 
equation 

Fk{s;xo,e,M) =det{Ak{s)). (5.11) 

Numerical results of computations of roots of the function Fk{s; xq, e, M) are shown on Figures 
1-5. Figure n show first few roots s of Fk{s] xq, e, M) with k = 1,3,5 versus xq for e = and 
M = 150. In agreement with Proposition |^ the roots converge as x ^ 1 to the values = c + n 
with o" = and n > 0. We can see that the convergence is excellent for k = 3 and k = 5 but it 
is worse for k = 1 in the sense that the roots at xq = 0.99 are still far from the values s„. This 
feature is explained by the decay of the eigenfunctions {^k, ^fc) of the system ()3.1|) on a; G [—1, 1]. 
Indeed, it follows from Proposition |21 that ~ (1 — x^)'^/^ and ^'a; ~ (1 — x2)i+°"/2 as x ^ ±1 
for e = and |A;| > 1. Therefore, the derivative of <I'fc(x) is bounded as x ^ ±1 for |A;| > 2 and 
unbounded for \k\ = 1. In the latter case, the power series expansions (|5.2() - H5.3|) diverge in the 
limit Xq ^ 1 and the numerical approximation is not accurate for xq close to 1. 

FigureEJshows first few roots s with k = 1,5 versus M for e = and xq = 0.9. We can see that the 
roots quickly converge to constant values, which are taken as approximations of real roots when 
M = 150 in the remainder of the figures. The numerical error for large values of M consists of 
three sources: truncation of the power series, root finding algorithms, and rounding entries of the 
matrix Ak{s) when a number xq with xq < 1 is evaluated at a large power Xq^. While the first two 
sources can be reduced to any desired degree, the last source represents an irremovable obstacle 
on getting accurate approximations when M gets large. 

Figure |21 shows the first six roots s versus k for e = 0, xq = 0.9, and M = 150. We observe 
two properties from this figure: the values of s becomes larger for larger values of k (e.g. the 
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Figure 1: First few roots s of Fk{s;xo,e^M) versus xq for e = and M = 150: k = 1 (circles), 
A; = 3 (stars) and A: = 5 (dots). 
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Figure 2: Convergence of roots s versus M for e = and xq = 0.9: k = 1 (circles) and = 5 
(dots). 
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eigenvalues // becomes more and more negative) and the roots s approach to the integer values for 
larger values of k even when xq = 0.9 is not close to xq = 1. 
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Figure 3: First six roots s versus for e = 0, xq = 0.9, and M = 150. 

Figure |1] show the first few roots s with k = 1,3 versus e for xq = 0.9 and M = 150. Although the 
roots are real for small values of e in agreement to Proposition El they coalesce for larger values of 
e. After two roots merge, they split into complex domain and complex values of s are not shown 
on Figure |^ It is seen from this figure that the roots with larger values of k coalesce for larger 
values of e. 

Figure [5] shows the spectrum of complex roots s with k = 1,3 for xq = 0.9, M = 150, and different 
values of < e < 12. The boundary of the stability domain ()5.1() is shown by the dotted curve. 
We can see that roots s remain in the stability domain after they bifurcate off the real axes. 



6 Numerical computations of eigenvalues for k = 

We approximate eigenvalues of the system (|4.1() with power series solutions explained in Section 5. 
The solution for ^o{x) and ^o{x) is represented by the power series (|5.2() - (|5.3|1 . where the starting 
coefficients (ao, &0) cq, do) are parameters, while the coefficients {am,bm, Cm, dm}m€N are defined 
uniquely from the recurrence equations. It follows from the ODE (j4.6j) that the set \^ciyyi, ^m}m$^N 
is uncoupled from the other coefficients but it is defined by the unknown value of the parameter s: 

_ {2m- s){2m + l + s)am-e{2m + l)bm , . 

"'"^^ - (2m + 2)(2m + l) ' ^^'^^ 

_ {2m + l- s){2m + 2 + s)bm-e{2m + 2)am+i „^ 
(2m + 3)(2m + 2) 
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Figure 4: First few roots s versus e for xq = 0.9 and M = 150: k = 1 (bolded curve) and k = 3 
(thin curve). 




Figure 5: Complex roots s for k = 1 (left) and k = 3 (right), xq = 0.9 and M = 150 when 
parameter e transverses in the interval < e < 12. The dotted curve shows the boundary of the 
stability domain (|5.1I) . 



23 



Given (ao, bo) and the value for s, the recurrence equation gives the value of ai and then the 
recurrence equation ()6.2() defines the value of bi, and so on. It follows from the first equation of 
the system ()4.1|) that the set {cm, dmjmeN is defined by the set {am, &m}meN but it is independent 
of s: 

_ (2m) (2m + l)cm + am ,„ „^ 

" (2m + 2)(2m + l) ' ^ ' 

{2m + l){2m + 2)dm + bm 
[2m + 3j(2m + 2) 

Finally, it follows from the second equation of the system (|4.1jl that there exist two initial equations: 

6o + eoo = -s{s + l)dQ, 
2ai + e6o = -2s(s + l)ci 

in addition to the system (|6.1|) - (|6.2|1 . When s / 0, we can solve the initial equations as 

6o = -eao - s(s + l)do, ci = ^, 

such that the only independent parameters are (ao,(io). We also note that the parameter cq is 
trivial since ^'0(2;) is defined up to the addition of an arbitrary constant. 

The boundary conditions in (|4.3|) lead to the equations: 

{^rn)cmx^^ = 0, (2"^ + ^)AraX^^ = 0. (6.5) 

m>0 m>0 

There exists a linear map from (oq, do) S parameterized by s G C to the sequence {am, bm, Cm, dm}m£N- 
Therefore, the boundary conditions ()6.5() are equivalent to the homogeneous system Aq{s)x = 0, 
where x = (ao,do)^ G and Aq{s) is a 2-by-2 matrix which depends on s G C, parameters xq 
and e, and integer M for truncation of power series. Eigenvalues ^ = —s{s + 1) of the system ()4.1() 
in (|4.3|) are equivalent to roots s of the determinant equation 

Fo(s;xo,e,M) =det(^o(s)). (6.6) 

Figure ini represents the first ten eigenvalues s versus xq for e = 1 and M = 100. In agreement with 
Proposition |SJ the roots converge to the integer values in the limit xq — > 1- Since the convergence 
of power series becomes slower with M for xq 7^ 1, there is a gap between the last numerical data 
and the value xq = 1. We also note that the numerical accuracy of the limiting eigenvalues (|4.16|) 
becomes worse for larger eigenvalues. 

Figure [7| represents the first ten eigenvalues s versus e for xq = 0.9 and AI = 100. It is obvious 
that the eigenvalues remain real in agreement with Proposition |3 

Figure |H1 represents the first seven eigenvalues s versus xq for e = 4 and two values of M = 100 
(dashed curves) and M = 1000 (solid curves). In agreement with Proposition llOl the roots converge 
to their limiting values which are not eigenvalues of the problem 1)4. 1|1 in space ()4.4() . We also note 
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limitations of the numerical methods based on truncations of the power series. True limits can 
only be recovered if too many terms of the power series are taken into accounts which leads to 
long computational time and large round-off errors of numerical computations. The effects of 
slow convergence and truncations of power series lead to coalescence of real eigenvalues and their 
splitting to the complex plane, which is not observed if the values of M are large enough. 
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Figure 8: Convergence of eigenvalues of the problem (|4.1|1 for e = 4 and two values of M = 100 
(dashed curve) and M = 1000 (solid curves). 



7 Discussions 



We have shown analytically that the stationary flow on the sphere is asymptotically stable whatever 
the Reynolds number may occur. This result is relevant for the flow of a viscous fluid (e.g. oil) 
over a sphere (e.g. a metal ball). We have also found that the linearized operator for symmetry- 
preserving perturbations has void spectrum in the energy space for sufficiently large Reynolds 
numbers. One can show by direct analysis that the full system (|1.5j) - (|1.7jl reduces to a scalar 
linear equation for symmetry-preserving (^-independent) solutions: 



dt sin ( 



(7.1) 



where Aq is given by (fm^ for A; = 0. When V0{e,t) = -%{e)e^\ the linear equation (TTT^ 
reduces to the linear eigenvalue problem (|2.18j) which has no eigenvalues in the space of square 
integrable functions Jq {"i!Q{9))^ sin 9 dO < oo when u < ^ {e> 2). Implications of this result to 
the well-posedness of the Cauchy problem for the linear time-dependent equation (|7.1|1 with < ^ 
remain unclear. 
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We have also shown analytically and numerically that the stationary flow on the truncated spherical 
layer is asymptotically stable and all isolated eigenvalues are real for small Reynolds numbers and 
complex for large Reynolds numbers. The eigenvalues are always real for symmetry-preserving 
perturbations. The truncated spherical layer can be used to model the ice melting in Arctics due 
to global warming, when the near-stationary flow of ocean water moves from Arctics to Antarctica. 
We note however that the model of two-dimensional Navier-Stokes equations on sphere considered 
in this paper does not include the Earth's rotation, the gravity force, and the location of continents. 
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